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Abstract: The spin period to precession period ratio of a non-axisymmetric spin-stabilized 
spacecraft, the Advanced Composition Explorer (ACE), was used to estimate the remaining mass 
and distribution of fuel within its propulsion system. This analysis was undertaken once 
telemetry suggested that two of the four fuel tanks had no propellant remaining, contrary to pre- 
launch expectations of the propulsion system performance. Numerical integration of possible 
fuel distributions was used to calculate moments of inertia for the spinning spacecraft. A Fast 
Fourier Transform (FFT) of output from a dynamics simulation was employed to relate 
calculated moments of inertia to spin and precession periods. The resulting modeled ratios were 
compared to the actual spin period to precession period ratio derived from the effect of post- 
maneuver nutation angle on sun sensor measurements. A Monte Carlo search was performed to 
tune free parameters using the observed spin period to precession period ratio over the life of the 
mission. This novel analysis of spin and precession periods indicates that at the time of launch, 
propellant was distributed unevenly between the two pairs of fuel tanks, with one pair having 
approximately 20% more propellant than the other pair. Furthermore, it indicates the pair of the 
tanks with less fuel expelled all of its propellant by 2014 and that approximately 46 kg of 
propellant remains in the other two tanks, an amount that closely matches the operational fuel 
accounting estimate. 
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1. Introduction 

ACE is a spin- stabilized spacecraft with two pairs of fuel tanks. In September 2013, one of four 
fuel tanks onboard ACE began to report increased temperatures. In May 2014, maneuvers began 
to periodically underperform by 10 - 20%. In November 2014, the other tank in the pair began 
reporting increased temperatures. In February 2015, maneuvers began to periodically 
underperform by as much as 50%. In this study, the ideal gas law is applied, showing the ratio of 
pressurant in one pair of tanks to that in the other pair is constant while fuel remains. A finite 
element model generates spacecraft moments of inertia over the time spanned by fuel remaining 
mass estimates. An FFT model using sun sensor data translates the history of moments of inertia 
into a history of predicted spin period to precession period ratios. A series of Monte Carlo 
searches were performed to show that the model can fit observed ratios but only if the initial fuel 
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distribution was uneven. The pressurant line design and unequal loading of pressurant hid this 
imbalance until anomalous propulsion system telemetry began to be received in 2013. All Monte 
Carlo trials that match (R 2 > 0.988) observed spin period to precession period ratios report one 
pair of tanks being empty in 2014. The Monte Carlo search performed in this study suggests that 
of the parameters explored, the actual spacecraft configuration deviates from the documented 
values only in terms of fuel distribution. The Monte Carlo search estimates 1.621 times as much 
propellant existed in the B tanks than in the A tanks at the time of launch. This corresponds to 
55.1% of the fuel being in the A tanks and 44.9% being in the B tanks at launch. 

2. Background 

2.1. Study Motivation 

The Advanced Composition Explorer (ACE) is a spin-stabilized spacecraft that was launched on 
August 25, 1997. ACE orbits the Sun-Earth LI libration point and provides near-real-time 
continuous coverage of solar wind parameters and solar energetic particle intensities. 

ACE performs attitude, orbit, and spin maneuvers using hydrazine thrusters. On average, ACE 
performs one maneuver per week. Four blow-down fuel tanks are onboard the spacecraft, split 
into two pairs (pair A and pair B). The propulsion system performed nominally for the first 16 
years of the mission. 

In September 2013 a small (0.01%) but unexplained increase in the spacecraft spin rate occurred. 
Also in September, the B1 tank began to report increased temperatures, as Fig. 1 illustrates. 



Figure 1. Daily average fuel tank temperature from telemetry 
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Figure 1 also shows the B1 tank temperature continued to increase to the point where in 
December 2013 the tank heaters — which had always been on for the length of the mission — 
began cycling on and off, causing the B1 tank temperature data to be less static. At this point, the 
mission began investigating a tank thermistor anomaly. 

However in May 2014, maneuvers began to periodically underperform by 10 - 20%. A 
retrospective investigation of the earlier increased temperatures and the response of tank B 1 to 
heater cycling as provided by Fig. 2 suggested tank B1 may be nearly devoid of fuel. This was 
initially met with some skepticism because the spacecraft propulsion system Critical Design 
Review (CDR) was interpreted at the time as implying the fuel should remain balanced onboard 
[ 1 ]. 
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Figure 2. Tank temperatures during heater cycling 

In November 2014, tank B2 began reporting increased temperatures and by December 2014 was 
reporting temperature data similar to tank Bl. The spacecraft design is such that all four tank 
heaters are either on or off. Thus, the temperature responses shown in Fig. 2 illustrate that both 
tanks in the B pair are acting in a manner consistent with having little to no fuel to act as a heat 
sink. Finally, in February 2015, maneuvers began to periodically underperform by as much as 
50%. 

Before any actions to mitigate the system performance were considered, this study was 
undertaken in November 2014 to determine if attitude data could be used to determine if the B 
tanks were nearly empty. This paper presents the results of that study. It also presents the results 
of a follow-on analysis which uses attitude data to determine the initial fuel distribution at the 
time of launch, along with other spacecraft mass properties. 

2.2. Spacecraft Layout 

The ACE propulsion system utilizes gaseous nitrogen to act as a pressurant in the blow-down 
hydrazine tanks. Fig. 3 is an illustration of the tanks, with the A-pair being the ones more closely 
aligned with the spacecraft body coordinate system (BCS) y axis. The system uses conispherical 
tanks, with the spacecraft spin forcing the fuel radially outward toward the tank nozzle. 
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Figure 3. Fuel tank configuration and BCS axes 

The fuel lines from all four tanks are interconnected. The pressurant lines, on the other hand, 
only connect tanks in a given pair, as is illustrated in the functional schematic provided in Fig. 4 
[2]. This arrangement allows one pair of tanks to have a different amount of pressurant than the 
other. A difference in the amount of pressurant obscures imbalances in fuel. While the four fuel 
tanks may report equal pressures, this does not equate to equal fuel masses; rather, it simply 
means more pressurant is present in the tanks with less fuel. Had both the fuel and pressurant 
lines been interconnected, the fuel would be equally divided if all pressures were equal. 



Figure 4. Fuel pressurant lines functional schematic 
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Once fuel is expelled from a pair of tanks and the corresponding fuel lines, pressurant escapes 
from the thrusters, causing marked underperformance. Since the fuel lines are interconnected, 
once a sufficient pressure difference exists between the two pairs of tanks, fuel from one pair can 
migrate to the empty fuel tanks and/or fuel lines. Once that fuel is expended, the cycle repeats 
itself. Hence, the design of the propulsion system allows for one pair of tanks to become empty, 
which then can cause periodic maneuver underperformance. This study uses attitude telemetry to 
determine if this is the cause of the telemetry and maneuver performance exhibited by ACE. 

Both star scanner and sun sensor telemetry is available from the ACE attitude control system. 
The sun sensor data is provided more frequently, though. The sun pulse time accuracy from the 
sun sensor is approximately 0.7 ms. Star scanner data is collected much more slowly. Normally 
the star scanner reports 4 stars every 64 s; for 20 minutes after each maneuver, 10 stars are 
reported every 16 s. Fine time resolution is needed in this study; consequently, sun sensor data is 
used and not star scanner data. Figure 5 provides diagrams showing the locations of the sun 
sensors, star scanner, thrusters. 
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Figure 5. Locations of sun sensors (SS), star scanner (ST), and thrusters 


3. Methodology 

In order to minimize precession induced by maneuver thruster pulses, engineers designed ACE 
to have a spin period to precession period ratio, P s:p , of 3:2, with a nominal spin period of 12 
seconds. In doing so, the perturbations induced by each thruster pulse act against the system 
precession and nutation [3]. Equation 1 relates the design principal moments of inertia to the spin 
period to precession period ratio. 

yyj Cpm 

?p recession y ^transverse, l x hrans ers 9,1 

The ratio P s:p can be determined using ACE sun sensor data after a maneuver. The induced 
precession of the spacecraft causes the sun sensor pulses to occur off of the nominal 12-second 
period. Plotting the variability of the time between sun pulses yields results like those shown in 
Fig. 6. Sun sensor data is used due to the finer time granularity of the data as compared to the 
star scanner. 
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Figure 6. Time between sun sensor sun pulses after maneuver on 2014/133 


The beat pattern shown in Fig. 6 occurs after each maneuver and is dependent on the spacecraft 
spin and precession periods. As Fig. 6 illustrates, the beats observed after the maneuvers have 
periods of approximately 6 s. This period, along with the spin period as measured by sun pulse 
timing, can be used to measure P s:p , using Eq. 2. 
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Figure 7 provides the history of this ratio from October 2001 to January 2015, calculated using 
Equation 2. When this ratio is plotted against estimated fuel remaining, a very strong linear trend 
is observed. The ratio becomes nearly constant when approximately 51 kg of fuel remain; this 
corresponds to the fuel estimate in August 2013, immediately before anomalous propulsion 
telemetry began to be reported. The estimated fuel is from a bookkeeping method, such as 
described in [4], and +3 kg error bars are provided, reflecting the uncertainty reported from the 
pre-launch blow-down curve. 
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A model of the spacecraft which determines the moments of inertia for a given fuel amount and 
distribution can produce an estimate of P s:p using Equation 1 . The results can then be compared 
to the observed ratio values obtained using Equation 2. When the modeled fuel distribution is 
very close to that actually onboard for a given fuel mass, the modeled and observed P s:p values 
should agree within some small tolerance. 

The value of the period ratio decreases only by 0.5% during the past 13 years, though. While 
Equation 1 provides an approximate relationship between the spacecraft moments of inertia and 
P s:p , simulation runs indicated the approximation was not accurate enough for usage in this study. 
Consequently, this paper introduces a new way to relate the moments of inertia to P s:p . 

First, a finite element model is used to calculate the moments of inertia of the distributed fuel 
mass; these are then added to the documented dry moments of inertia to determine the spacecraft 
moments of inertia for a given fuel mass and distribution. Next, a numeric model of Euler’s 
equation and a fast Fourier transform (FFT) are used to generate a look-up table of P s:p values for 
given triplets of principal moments of inertia. The ideal gas law is used to determine the fuel 
distribution during the course of the mission for a given pre-launch fuel distribution. Finally, a 
Monte Carlo analysis is performed to determine what initial distribution of fuel generates period 
ratios that most closely resemble observed data, taking into account uncertainty in quantities 
such as the dry moments of inertia and tank nozzle locations. 

The next sections of this paper provide additional model details. But the overarching goals of this 
study are to use these models to determine the fuel balance onboard and to determine if the B 
tanks became empty in late 2014. A secondary goal of this study is to use a Monte Carlo search 
to estimate the dry moments of inertia and the tank nozzle locations more accurately than 
documented values available to the mission. 
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4. Model Specifics 

4.1. Fuel Tank Finite Element Model 

The fuel mass in each tank is modeled using a 60 x 60 x 60-element grid, with each element 
enclosing 0.21 cm 3 of volume at the start of a nominal run. The location and geometry of the fuel 
tanks are provided by mechanical drawings available in the ACE propulsion subsystem critical 
design review (CDR). The drawings are not complete and visual inspection was used to estimate 
the fuel tank nozzle locations. The initial fuel mass and the dry mass properties are taken from 
the ACE pre-ship review [5], Numeric integration of the finite element model is used to 
determine the fuel — and hence the overall — moments of inertia. 

Propellant residing in the fuel lines is neglected. Slosh is also neglected, along with any other 
phenomena that would deform the propellant/pressurant interface. Hence, the fuel in each tank is 
modelled as a slug pinned against the outermost surface of the conispherical tank with the 
propellant/pressurant interface being parallel to the spacecraft spin axis. Additionally, the 
hydrazine density is treated as a constant value of 1.08 g/cm 3 . 

Due to the shared fuel and pressurant lines between the two tanks in a pair, each tank in a pair is 
treated identically in terms of the amount of pressurant and propellant present. In subsequent 
equations, each pair of tanks is treated as one unit. 

Until 2013, pressure readings from the A tanks and B tanks matched within 1 least significant 
bit. Temperature readings also agreed well, and even after 2013, temperature differences were at 
most 5 K. Consequently, when modelling the nitrogen pressurant in the A and B tanks, the 
pressures and temperatures are treated as being equal. Since nitrogen behaves as an ideal gas at 
the temperatures and pressures observed in the tanks, the ideal gas law can be applied and 
simplifies to that shown in Equation 3, where n is the number of moles of nitrogen pressurant 
and V is the volume of the pressurant. 


n B tanks ^StQnfcS-pr&ssurant 

' (3) 

n A tan ks v Ata n ks_p ressu ran t 

Since the pressurant is in a closed system — at least until nitrogen begins to escape due to little 
fuel remaining — the number of moles of nitrogen in the A and B tanks remains constant. Thus, 
the ratio of the nitrogen volume in the B tanks to the A tanks at the time of launch remains 
essentially constant through the course of the mission. Equation 4 shows the definition of a 
variable nitrogen volume ratio, B to A, (Vn2_bwa). 

V N2 BtoA — — constant ( 4 ) 

^Atan ks_p rsssu rant 

Using the fuel mass and the density of hydrazine, the total hydrazine volume can be easily 
computed. The remaining tank volume, by definition, must be occupied by the nitrogen 
pressurant. The ratio VN 2 _BtoA is then used to determine the volume of pressurant in the A tanks 
and in the B tanks. Equation 5 describes how to calculate the volume of the pressurant in the A 


8 


tanks; the volume of pressurant in the B tanks can be determined by applying the ratio VN 2 _BtoA. 
The distribution of the fuel between the A tanks and B tanks follows directly from this, for the 
tank volume not consisting of pressurant consists of fuel. 


v A 


Atanks jpr&ssu ran t 


m hydrazin 

(Ttcmfcs rotaf - n J 

~ ^hydrazine 

( 1 + V N 2 BtoA ) 


( 5 ) 


The finite element model is used in the following fashion to determine spacecraft moments of 
inertia. A test value of Vn2_bioa is selected randomly. For each estimated remaining fuel mass in 
Fig. 7, the distribution of the fuel between the A tanks and B tanks is determined. The 
assumption that individual tanks in a pair are identical is then used to determine the fuel mass in 
each of the four tanks. Finally, the finite element model is used to determine the moments of 
inertia induced by the fuel. This is added to the dry moments of inertia to determine the total 
spacecraft moments of inertia for the given fuel mass and VN 2 _BtoA. 


4.2. Fast Fourier transform model for P s:p 


Euler’s equations for torque-free motion provide a means to translate from spacecraft moments 
of inertia to measurements of the spin and precession periods. Equations 6-8 list these 
equations, without any simplifications being applied (treated as asymmetric and precession and 
nutation are not neglected). One assumption being used, though, is that the net moment due to 
fuel slosh is zero, as are any external torques such as that due to solar radiation pressure. 

/jCPj (/ 2 — ^3)^2 ^3 (6) 

I 2 d>2 = (Jl — h ) ^3 ( 7 ) 

^3^3 (Jl ^2)^1 ^2 (8) 


Equations 6-8 were used to numerically propagate a model of ACE, using an initial spin rate. 
Principal moments of inertia were selected combinatorially from a 1-kgm 2 grid of values that 
spanned all moments of inertia observed in this study. Propagations used a Runge-Kutta method 
(Dormand-Prince 8(5, 3)), and an initial spin period of 12 s. Without loss of generality, here C 03 is 
the spin rate about the BCS Z axis and 001 and 002 are rotations about the BCS X and Y axis, 
respectively. An initial angular velocity about coi similar to that seen immediately after 
maneuvers is also used. All other angular velocity and accelerations are 0 initially. 

In this method, data relating the location of the spacecraft body axes to the location of inertial 
axes are used indirectly to measure spacecraft spin and precession. Figure 8 illustrates the layout 
of these axes. The times when a body transverse axis crosses the inertial x-y plane are recorded; 
this provides a measure of the rate of precession. Similarly, the times when that same transverse 
axis crosses the inertial x-z plane are recorded to measure the rate of spin. 
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Figure 8. Simplified diagram of body and inertial frames used in FFT model 

A fast Fourier transform was performed on these times to determine the periods of the resulting 
spin and precession. Long propagations (8,192 s = 2 13 s) were used to obtain high-precision 
frequency measurements from the FFT; the resulting measurement precision of P s:p is 0.0016. 
Figures 9 & 10 provide sample output, showing the strong frequency responses that correspond 
to the spin and precession frequencies (from which the periods follow). 
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Figure 9. Sample FFT frequency response for spin (f = 0.08350 Hz, T = 11.98 s) 
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Figure 10. Sample FFT frequency response for precession (f = 0.1125 Hz, T = 7.968 s), also 
showing response for nutation (f = 0.04102 Hz, T = 23.38 s) 

In order to decrease simulation run time, the FFT model was initially run repeatedly with the 
aforementioned sample space of I-triplets to generate a lookup table relating moments of inertia 
to P s:p . The lookup table uses a grid spacing of 1.0 kg-m 2 for Ii, h, and I 3 . Once the finite element 
model generates moments of inertia, the lookup table is used to supply the corresponding P s:p , 
rather than having to run a time-consuming FFT each time. This is especially important for the 
Monte Carlo analysis performed in this study. 

4.3. Model Summary 

Figure 1 1 supplies an overview of the complete ACE model used in this study. First, a value of 
V N2_BtoA at the time of launch is selected, which is used to determine the fuel distribution over 
time. Estimated fuel masses over the past 13 years are then collected. Both of these are supplied 
to the finite element model, which generates spacecraft moments of inertia over the time spanned 
by the fuel mass estimates. The inertia matrices are diagonalized to determine the history of the 
principal moments of inertia. Next, the FFT model translates the history of moments of inertia 
into a history of predicted P s:p values. Finally, the predicted ratios are compared to the observed 
ratios and the coefficient of determination, R 2 , is calculated to determine the degree of fit. 
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Figure 11. ACE model overview 

Two Monte Carlo analyses are performed in this study. First, the value of VN 2 _BtoA is allowed to 
vary over a wide range and uncertainty values are assigned to dry moments of inertia and the 
tank locations. A wide range of VN 2 _BtoA is initially used since the fuel distribution is unknown. 
The Monte Carlo results are used to determine the likelihood that the initial distribution of fuel — 
and by Equation 4, the fuel distribution for the entire mission — was asymmetric. Specifically, 
trials using values of VN 2 _BtoA that generate predicted P s:p values that closely match observed P s:p 
values are collected. From this, the likelihood of the B tanks being currently empty is assessed by 
examining the fuel distribution in 2014 modeled by the collected successful trials. 

Next, after identifying the range of VN 2 _BtoA that is consistent with observations, a second, finer 
Monte Carlo search is performed. This search is executed to identify sets of values for VN 2 _BtoA, 
dry moments of inertia, and tank locations that generate P s:p values that most closely match 
observations. 

4.4. Selection of Monte Carlo perturbations 

The ACE CDR states that the launch configuration moments of inertia have an accuracy of 1.5%, 
without stating whether these are 1 -sigma values, or 95% confidence intervals, or any other 
measure. For the purposes of this study, all dry moments of inertia are allowed to vary by +/- 
4.5%. The locations of the fuel tank nozzles were identified using the mechanical drawings 
presented in the ACE propulsion portion of the ACE CDR. However, these drawings provided 
only estimates of the nozzle locations within Vi inch of a nominal location. In order to account 
for uncertainty in this estimate, and uncertainty in the difference between the design placement 
and the actual placement, the tank nozzle coordinates are allowed to vary by +/- 1 inch. Lastly, 
those tanks nozzle coordinates were documented in the propulsion coordinate system (PCS). 
This system is aligned with the BCS, but the origin is offset (ZBCSoffset) in the +Z direction from 
the BCS. In order to write the tank coordinates in the BCS frame, ZBCSoffset was first determined. 
Hardware for which coordinates could be found in both the BCS and PCS frames were used to 
determine this offset. Again, though, to allow for uncertainty in actual hardware placement — and 
to be consistent with the uncertainties mentioned thus far — ZBCSoffset is also allowed to vary by 
+/- 1 inch. 
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Since the first Monte Carlo analysis is being performed as a search for possible values of 
V N 2 _BtoA s rather than determining the sensitivity of VN 2 _BtoA to modeled input parameters, all 
uncertain parameters were sampled from uniform distributions with the aforementioned limits. 
V N2_BtoA was also sampled from a uniform distribution, spanning 0.5 - 2.5. Table 1 summarizes 
the nominal values and uncertainty bounds for the perturbed variables used in the first Monte 
Carlo analysis. 


Table 1. Sample Space used in First Monte Carlo Analysis 


Quantity 

Nominal Value 

Sample Space 

V N2_BtoA 

1.0 

U[0.5,2.5] 

Ixx,dry 

324.44 kgm 2 

U[nominal - 4.5%, nominal + 4.5%] 

Iyy,dry 

249.59 kgm 2 

U[nominal - 4.5%, nominal + 4.5%] 

Izz,dry 

436.16 kgm 2 

U[nominal - 4.5%, nominal + 4.5%] 

Ixy,dry 

-0.14 kgm 2 

U[nominal - 4.5%, nominal + 4.5%] 

Ixz,dry 

0.03 kgm 2 

U[nominal - 4.5%, nominal + 4.5%] 

Iyz,dry 

0.01 kgm 2 

U[nominal - 4.5%, nominal + 4.5%] 

ZBCSoffset 

0.280 m 

U[nominal - 1”, nominal + 1”] 

A1 nozzle location 

(-0.569, -0.236, 0.109) m 

x,y,z e U[nominal - 1”, nominal + 1”] 

A2 nozzle location 

(0.569, 0.236, 0.109) m 

x,y,z e U[nominal - 1”, nominal + 1”] 

B 1 nozzle location 

(0.236, -0.569, 0.109) m 

x,y,z e U[nominal - 1”, nominal + 1”] 

B2 nozzle location 

(-0.236, 0.569, 0.109) m 

x,y,z e U[nominal - 1”, nominal + 1”] 


5. Results 

5.1. Initial Monte Carlo search results: Are the B tanks empty? 

For each Monte Carlo trial, P s:p was calculated for each estimated fuel mass used in the 
observational data. To compare each trial’s history of ratios to the observed history of ratios over 
the past 13 years, the coefficient of determination, R 2 , was calculated. 

Initial Monte Carlo runs were executed and trials with large values of R 2 (R 2 > 0.9) were 
examined. Trials with R 2 > 0.988 were judged to have very close agreement between predicted 
and observed ratios. This lower bound is used to define when a trial is a ‘success.’ A value of 1 is 
assigned to each successful trial, and 0 is assigned to all others. 

With the limit now selected, a 12,000,000-trial Monte Carlo was performed and 175 successful 
trials were identified. For all successful trials, VN 2 _BtoA is contained in [1.531, 1.706]. To 
indirectly measure the likelihood that some other region of the sample space would have yielded 
successful trials, a confidence interval was calculated. The actual, population ratio of the space 
that yields a success to the entire sample space is p. The Monte Carlo analysis provides a sample 
estimate, pn, of 1.46xl0" 5 . The Clopper-Pearson 95 th percentile confidence limits for this 
binomial distribution is [1.25 xlO" 5 , 1.69 xlO" 5 ]. Given the relative tightness of this range 
compared to the sample estimate, it is unlikely that other regions in the sample space would yield 
successful fits to the observed data. 
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To visualize the level of agreement obtained by the successful trials, Fig. 12 provides a plot of 
the predicted and observed P s:p for trial 1709670, for which R 2 = 0.998. 



• Measured Ratios 

• Modeled Ratios 


92 82 72 62 52 42 

Estimated Fuel Remaining (kg) 

Figure 11 . Modeled and measured F s:p lor trial 170^070, K" = U.yys 


The range of VN 2 _BtoA for the successful trials indicates the fuel was not evenly distributed among 
the tanks at the time of launch. Lending further credence to this finding, the ‘elbow’ present in 
the observed and simulated data for the successful trials corresponds to approximately one month 
before the B1 tank began reporting increased temperatures. Hence, the elbow in the data occurs 
at nearly the same time the B1 tank was believed to have become nearly devoid of fuel. 


By contrast, trials in which the value of VN 2 _BtoA was nearly 1.0 (i.e., the fuel is divided equally 
in the tanks and fuel currently remains in all tanks) do not perform well, particularly with respect 
to fitting observed data once B tank telemetry became anomalous. Figure 13 provides output 
from the trial with VN 2 _BtoA ~ 1 that reported the highest value of R 2 (VN 2 _BtoA =1.002, R 2 = 
0.978). The uncertainty values of the measured ratios originally shown in Fig. 7 are repeated 
here. The key feature in the observed history of P s:p — the flattening of the ratio when 
approximately 51 kg of fuel remains — is missing from this trial and all other trials with VN 2 _BtoA 
~ 1. While the uncertainties in the measured data allow for uncertainty in the slope of the linear 
trend of P s:p while more than 51 kg of fuel remains, they do not allow for that trend to continue 
past the elbow point. Thus, evenly distributed fuel at the time of launch is incompatible with 
observed data. 
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• Measured Ratios 

• Modeled Ratios 


Figure 137 Modeled F s:p tor best fitting trial with VN2_BtoA ~ 1 (K^ = U.y78, VN2_BtoA = 1.002) 
as compared to measured history of P s: p with uncertainty bars 


5.2. Refined Monte Carlo search results: Identifying actual spacecraft configuration 

During the second Monte Carlo analysis, Vn2_bioa was sampled from a uniform distribution 
spanning the range of successful values identified in the first analysis, [1.531, 1.706]. In the set 
of successful trials from the first Monte Carlo analysis, minimal correlations were present 
between all tank location input variables. Also, for those successful trials the coordinate values 
essentially spanned their sampled range. Hence, for the second Monte Carlo analysis, all input 
values related to tank locations were sampled in the same fashion as during the first analysis. 


In the set of successful trials from the first Monte Carlo analysis, values of I xx ,dr y and I yy ,dry also 
spanned their sample space and showed negligible correlation to each other (r = -0.056). From 
that set of successful trials I zz ,d ry did not span the sample space (instead spanned 416 - 448 kgm 2 ) 
and visually did not appear to be uniform in distribution. Tests using the sample second, third, 
and fourth moments were not able to identify the distribution of I zz> d ry , though [6]. Although I xx ,d ry 
and I yy ,dr y exhibited strong correlation with I zz ,d ry (r = 0.869 and 0.445, respectively), the lack of a 
distribution to use for I zz ,d ry is an impediment to using the correlation between the dry moments 
of inertia when sampling Ixx,dry? Iyy,dry, and Izz,dry for the second Monte Carlo analysis. 
Consequently, I xx ,d ry and I yy ,dr y were sampled uniformly in the same fashion as the first analysis 
and I zz ,diy was sampled uniformly from [416, 448] kgm 2 . 

The second analysis used 4,000,000 trials. The value of R 2 is again used to measure the degree of 
fit between the predicted P s:p and observed P s:p values. The 0.01% of the trials with the highest 
value of R 2 were collected and represent the best fit set. Table 2 presents the median, x, and 
sample standard deviation, s, for the input variables from the best fit set, along with the nominal 
values. Figure 14 provides the output from the trial with the highest value of R 2 (R 2 = 0.998). 


15 



Table 2. Best Fit Set Estimates of Actual Spacecraft Configuration, Compared To Nominal 


Val 

ues 

Quantity 

Best fit set x 

Best fit set s 

Nominal value 

V N2_BtoA 

1.621 

0.027 

1.0 

Ixx,dry 

327.89 kgm 2 

7.70 kgm 2 

324.44 kgm 2 

i— i 

a. 

246.47 kgm 2 

6.32 kgm 2 

249.59 kgm 2 

Izz,dry 

43 1 .04 kgm 2 

7.58 kgm 2 

436.16 kgm 2 

ZBCSoffset 

0.292 m 

0.011 m 

0.280 m 

A1 nozzle location 

(-0.576, -0.238, 0.113) m 

Sx = Sy — Sz — 0.014 m 

(-0.569, -0.236, 0.109) m 

A2 nozzle location 

(0.576, 0.236, 0.116) m 

Sx=s y =s z = 0.014 m 

(0.569, 0.236, 0.109) m 

B 1 nozzle location 

(0.239, -0.567, 0.108) m 

Sx = Sy — Sz — 0.014 m 

(0.236, -0.569, 0.109) m 

B2 nozzle location 

(-0.236, 0.569, 0.110) m 

Sx = Sy — Sz — 0.015 m 

(-0.236, 0.569, 0.109) m 



• Measured Ratios 

• Modeled Ratios 


92 82 72 62 52 42 

Estimated Fuel Remaining (kg) 

Figure 14. Modeled and observed F S : P for trial 7S7U22, = U.WS 


All nominal values are within one sample standard deviation of the best fit set median values, 
except for VN 2 _BtoA, which is restricted to a tight range that does not include 1.0 (x= 1.621, s = 
0.027). Consequently, results from this study’s model are consistent with the documented 
spacecraft configuration with one exception. The pressurant — and hence the propellant — was 
unevenly distributed onboard at launch. The median value of 1.621 for VN 2 _BtoA corresponds to a 
fuel distribution of 55.1%/44.9% between the A and B tanks, respectively, at the time of launch. 


6. Conclusion 


Based on the results of this study, it is very likely that at the time of launch the fuel was not 
evenly distributed among the four fuel tanks. The pressurant line design and unequal loading of 
pressurant hid this discrepancy until anomalous propulsion system telemetry began to be 
received in 2013. All Monte Carlo trials that match (R 2 > 0.988) observed spin period to 
precession period ratios report the B tanks as being empty in 2014. All values perturbed in the 
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Monte Carlo searches performed in this study estimate spacecraft configuration values that are 
within one sample standard deviation of pre-launch documented values with one exception — the 
fuel distribution. The Monte Carlo search estimates 1.621 times as much propellant existed in the 
B tanks than in the A tanks at the time of launch. This corresponds to 55.1% of the fuel being in 
the A tanks and 44.9% being in the B tanks at launch. Based on the findings of this study, the 
ACE mission closed the B-side latch valves on March 17, 2015 to prevent pressurant from 
escaping through the thrusters which thereby causes maneuver underperformance. Since doing 
so, maneuvers have returned to historically nominal performance, with only one maneuver (one 
with a magnitude much lower than typical maneuvers) underperforming by more than 5%. 
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